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Abstract 

Different aspects of the Gribov-Glauber approach for calculation of inclusive pion double 
charge exchange (DCX) on nuclei are investigated. Recently we have shown that inelastic 
rescatterings (IR) with two (and more) pions in the intermediate states give an important 
contribution to the process of DCX at energies above ~ 0.6 GeV . In this paper we use 
the one pion exchange model to study in details amplitudes of two pion production. This 
allows us to verify theoretical assumptions made in the previous paper and to predict 
cross section for the forward inclusive pion DCX at energies > 1 GeV, where IR dominate 
over the conventional DCX mechanism of two sequential single charge exchanges. 
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1 Introduction 

The inclusive pion double charge exchange (DCX) reaction on nuclei, e.g. 

tC + A(Z, N) -> 7T+ + X (1) 

as a two-step transition on two like nucleons (protons) is described by the 
diagrams of three types (Fig. 1(a), (b), (c)) which correspond to the Glauber 
picture of elastic (a), quasielastic (b) and inelastic (c) rescatterings. 

The conventional DCX mechanism of two sequential single charge ex- 
changes (SSCX), i.e. the elastic rescattering, is found to describe [1], in gen- 
eral, experimental data on pion DCX at energies of the meson factories up 
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to 0.5 GeV . It predicts [2] a strong decrease of forward angle exclusive pion 
DCX at incident kinetic energies Tq > 0.6 GeV due to the rapid energy drop 
of 7r~ — > 7T° transition rate. The SSCX mechanism leads [3] to analogous be- 
haviour of inclusive forward DCX rate too. 

Inclusive reaction (1) is a unique process where, as it was shown in [4], 
the effect of the Glauber inelastic rescatterings (IR) becomes very important 
already at GeV energies. This effect can be expected from a comparison of 
experimental total cross sections for two competing processes, 7T~p — > 7r°n 
and 7T~p — > 7T + 7r~n (Fig. 2, where data are given from the compilation [5]). 
The study of the forward angle reaction (1) by ITEP experiments [6] on 16 
and 6 Li in a kinematics which forbids real production of additional pion has 
demonstrated a relatively weak energy dependence of its rate and considerable 
excess over SSCX mechanism prediction [6] around Tq = 1 GeV. So the con- 
ventional elastic SSCX picture of inclusive pion DCX has to be modified by 
an important new mechanism of IR at GeV energies. 

Our goal is to elaborate recently suggested Gribov-Glauber approach to 
inclusive pion DCX reaction [4] on 16 0. Below we present an updated analysis 
of some theoretical aspects of our model and clarify simplifying assumptions 
we have used. Then we estimate quantitatively upper and lower limits for in- 
clusive pion DCX cross section in the energy region of Tq = 1 GeV to 4 GeV 
to be tested in future experiments. Our working scheme for the calculation of 
7r — > 2jt amplitude is the one pion exchange (OPE) model which we check by 
a comparison of theoretical cross section with experiment as well as distribu- 
tions on two-pion effective mass for the process n~p — > 7r + 7r~n. We compare 
our results with predictions of other models and in particular of the meson 
exchange currents (MEC) model [7]. 

2 Gribov-Glauber picture of inclusive pion DCX 

In this Section we shall discuss some aspects of the formalism for a calcu- 
lation of the amplitude of the DCX process n~pp — > 7T + nn (Fig. 3). Following 
Gribov [8] it is convenient to introduce an integration over Mjy, square of mass 
of the intermediate state in Fig.l. It is related to integration over Fermi motion 
of initial protons and momenta of final neutrons. These momenta are limited 
by experimental conditions. For example in the ITEP experiment on 6 Li and 
16 [6] the inclusive DCX reaction A(tt~ ,tt + )X was measured at kinetic en- 
ergies Tq = 0.6, 0.75 and 1.1 GeV ((6) « 5°) with kinematical condition 
AT = Tq — T < m n (T is the kinetic energy of the outgoing pion) in order to 
exclude additional pion production. 
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The kinematical region indicated above strongly limits momentum trans- 
fers t\{2) to neutrons in the transitions TT~p —> Hn and Hp — > tt + u (see 
Fig. 3). Mh can be very large at high incident energy, Eq, and increases as 
M]j ~ IijinEq. The value of Mjy varies also due to the Fermi motion. It 
is easy to obtain an estimate of this variation taking into account that at 
high energies, E > m N , \t^\ « (M| - mj) /4£ 2 < ^ 2 < 2m A rAT 1 ( 2 ) and 

(M 2 H ) max « 2Eo(2AT l m iV ) 1/2 . For energy £ - 1 GeF, {M 2 H ) max is equal to 
0.5 GeV 2 . As it was explained in Ref. [4], the integral over can be written 
either over the real axis or (due to analyticity in M|) as a sum of the ab- 
sorptive part and integral over semicircle at \M^\ ~ (Mj I ) max . At high enough 
energies i?o the last contribution can be neglected if the corresponding ampli- 
tude decreases faster than 1/Mjj at large Mjj. We gave arguments [4] that the 
amplitude of Fig. 3 should satisfy this property. In this case the DCX amplitude 
is proportional to the discontinuity in Mjj of the amplitude in Fig. 3, and using 
unitarity we obtain the representation equivalent to the diagrams of Fig.l. 

Contribution of a single particle intermediate state (tt°, 77°) or narrow 
resonance (u, ...) to the DCX cross section can be written in the form 

Let us comment on contribution to DCX of different diagrams of Fig.l 
that was taken in Ref. [4] additively. The diagram (6) gives a small correction 
to the Glauber elastic rescattering of Fig. 1(a) due to the fact that 7r — — ^ 77 
amplitude is relatively small (see Fig. 2) while the production of two pions is 
a very important competing mechanism, especially for Tq > 0.6 GeV. As for 
diagrams (a) and (c), if tt — > 2tt amplitude is dominated by the pion exchange 
(see Fig. 4), then their interference is absent. Really, the invariant form of 
the nucleon vertex in Fig. 4 is u^u, while the amplitude of the transition 
7r~p — > 7T°n is proportional to u(A + Bq)u, where q = p n - + p n o. So the 
interference term is Tr(z3i + m^lhip'i + itin)(A + Bq) = 0. 

Note that for the s-wave tttt production which is a dominant inelastic 
process at energies Eq < 1 GeV the interference is absent by the same reason. 
For the s wave (and in general for all even waves) the same expression (2) is 
valid. However for realistic case of tttt production, when both even and odd 
orbital momenta are important, the amplitude of tt~ — > tt + transition (contrary 
to the diagonal tt~ — > tt~ transition) can not be written [4] in the form (2). 
This is especially clear in the tt exchange (OPE) model (Fig. 4) where DCX 
amplitude is expressed in terms of the tt~tt + — > tt + tt~ amplitude (or backward 
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elastic 7r~7r + scattering amplitude). 

In this paper we shall use the OPE model in order to study the problems 
outlined above: is it possible to neglect the contribution of the large semicircle 
at energies 1 GeV and what is the difference between integrals of forward and 
backward 7r~7r + scattering amplitudes? In this way we shall find corrections 
Tjj to simple expressions for DCX cross sections 




where 



,exp ,exp \ 2 

C lmax i2 ( ]\/f 4. \ I r lmax » 

V h min (M) aMal1 ' J 



(3) 



(4) 



r^o = (cr^t^o/ a^i^-) 2 • r n+n - (5) 

introduced in Ref.[4] and shall obtain lower and upper bounds for these cross 
sections at higher energies. 

3 Inelastic rescatterings and comparison with ITEP experiment 

In this Section we shall calculate the contribution of IR to DCX using 
Eqs.(3)-(5) and experimental data on particle production in reactions 

Ti p — > Hfl, H = TT°, TT+7T-, TtV, ... (6) 

The quantities T# in Eqs.(3)-(5) were determined using experimental data 
for the following hadronic states H: tt° [9], 7r + 7r~ and 7r°7r° [10]. While inte- 
grating over t and M 2 , experimental constraints on AT are taken into account 
(see details in [4]): \t m \ < 2m N {AT max /2) Q. 

In the Table 1 we show the updated summary of energy dependence for 
forward inclusive DCX reaction cross sections of tt~ on 16 obtained in exper- 
iments at ITEP [6] and calculated via Eqs.(3)-(5). The first line in the Table 1 
contains the values of da ll DCX / d^l calculated [3,6] in the framework of the elas- 
tic SSCX mechanism without taking into account effect of the Fermi motion 
(FM) on magnitude of the 7riV-interaction energy in the nucleus. In the values 
d&DCX I d£l of the second line the FM effect Q is included. The resulting values 

1 J In principle it is necessary to integrate over t\ and £2, but taking into consideration some uncertainty of 
the theoretical estimates as well as of the experimental information available we limited ourselves with this 
simple relation. 

2 -* We thank L.Alvarez-Ruso and M.J.Vicente Vacas for preparing the version of the code [3] which takes into 
consideration the Fermi motion. 
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are daocx/dVL = da^ cx /dQ + (IY^- + T n o w o) da^ cx /dQ are given at the 
fourth line. Due to a sharp change of DCX rate at energies of interest the FM 
effect gives smearing of the dip-bump structure, that in average leads in some 
increase of the values of DCX cross section. Note, that all values in the Table 
1 correspond to the outgoing momentum spectra integrated over the region of 
AT = to 140 MeV. 



4 Inclusive pion DCX in the energy range of 1 — 4 GeV 

We see from the experimental data compilation of Fig. 2 that the cross 
section values of the reaction (6) with H = 7t + tt~, 

7T~P — > 7T + 7T~n, (6') 

exceed ones with H ^ 7r + 7r~ not only at 1 GeV but also at higher energies. 
So it is reasonable to assume that IR with two pion intermediate states will 
dominate inclusive pion DCX at least up to 4 GeV . In this case the diagram 
of Fig. 4 can be used to obtain a quantitative estimate of DCX rate within 
the OPE model. Note here that, besides of the well-known good description of 
experiment for (6') reaction, the OPE model is very convenient to our goals 
because of its factorization feature. This permits to separate M 2 and t depen- 
dences of the (6') reaction amplitude, A n ^2n(M 2 , t), and to take into account 
limitations on their values in the experiment [6] . 
We start from the following formula 



da™ 



DCX oc 



dn 



( M 2 (t max ) t max \ 2 

/ dM 2 J dtlmAiZT d (M 2 ,t) (7) 

V 4m2 t min (M 2 ) J 

equivalent to Eqs.(3)-(5). Then we use the OPE model [11] to calculate a 
difference between imaginary parts of forward and backward amplitudes and 
to estimate the real part of the amplitude. 

In OPE model it is possible to express the amplitude of the process (6 ; ) 
in terms of the tttt scattering amplitude, 

A^ 2n (M 2 ,t) = F 2 (t) A nn (M 2 ) (8) 

where 

F 2 (t) = 2G 2 [\t\/(t - ml) 2 } exp[2R 2 (t - ml)}, R 2 = 1.92 GeV~ 2 , (9) 



5 



G is the 7riViV-coupling constant. 
(a) OPE model for DCX 

As it was shown in [11] the OPE model calculations are in a good agree- 
ment with the experimental data on the reaction (6') above 2 GeV. So for 
calculation of the cross section for this reaction we shall follow Ref.fll], namely: 

where G 2 /4tt =14.6, Q(s,ml,m^) = ^s 2 - 2s(m\ + mf) + ( m i ~ ^i)/ 2 \/s, 
M is the mass of two-pion state, 

^ +7r -(M 2 ) = 2?r f | | 2 efe, 2 = cos#, (11) 

and the 7T7T scattering amplitude is equal to 

MM\z) = ^L^Z(2l + l)ll + (-iy+']f[(M')P,(z). (12) 

By virtue of the unitarity condition we have: 

Im^(M 2 , z = 1) = 2Q(M 2 , m 2 , m 2 )M^(M 2 ). (13) 

The partial- wave amplitudes fi (M 2 , z) for orbital angular momentum / and 
isospin / are expressed via the phase shifts, Sf , and elasticities, rjf: 

//(M 2 ) = I[^)e^) - 1], = !/> + 1/0 + 1/1. (14) 

It should be noted that the equation (7) takes into account both tt + 7t~ 
and 7r°7r° intermediate states. 

In the following we will use the complete set of amplitudes (phase shifts 
and elasticities) with / = 0, 2 for / = 0, 2 and with I = 1, 3 for / = 1 in the 
range of M up to «1.8 GeV from Ref.[12] where fixed-/; and fixed-w analyticity 
in conjunction with energy independent phase-shift analysis gave the single 
solution. The following features of the model confirm that its predictions are 
reasonable. 

(1) The M 2 dependence of the total cross section of 7r + 7r~scattering calculated 
according to Eq.(ll), (13), (14) is in a good agreement with the data from 
Ref.[13] (see Fig.5). 

(2) The value of the parameter R 2 = 1.918 GeV~ 2 that we use was obtained 
on the base of our analysis of the t dependence of the yield for the reaction (6') 
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averaged over the energy region from 2 to 3 GeV taken from [11]. The total 
cross section of this reaction calculated via the integration of Eq.(lO), 



in Fig. 6. It is seen that the model used (full circles) reproduces experimental 
energy behaviour of a^-p^+^-n only above «2.5 GeV. Some underestimate of 
the absolute value of the cross section decreases as To increases. 

In Fig. 6 the cross section of the reaction (6'), obtained with constraints 
of tmax = tmax = 0.135 (GeV/c) 2 corresponding to the experimental limita- 
tion AT < 140 MeV, is also presented. The experimental values (crosses and 
empty star) were calculated from the data of Ref.[10a] on M 2 and angu- 
lar cross-section dependencies of the reaction (6') (crosses) and taken from 
Ref.[14] for the interval 1.5 < tfm\ < 8 (empty star). The OPE model used 
(triangles in Fig. 6) agrees with experiment starting already from 1 GeV. Note 
that below Tq = 1 GeV the OPE model predictions strongly deviate from ex- 
perimental data. This is probably related to a substantial role of s-channel 
resonance production in this region. 

(3) The experimental distributions on the two-pion mass (M or M 2 ) for inci- 
dent momentum values of 1.343 [10], 1.59 [14], 2.26 [15], and 4 GeV/c [16] are 
compared with the calculated ones (see Figs. 7-10). It is seen from Fig. 7 that 
the model predicts too small cross section in the low M 2 region at the lowest 
energy. On the other hand, at higher energies the OPE model reproduces ex- 
periment rather well (including p and /2(1270) resonances as it is seen in Figs. 
9, 10). For our purpose it is essential to have a good description of experiment 
in t interval close to t < 0.135 (GeV/c) 2 . It follows from Fig. 8 that for the 
region of 1.5 < t/m 2 < 8 an agreement with experiment is better than for the 
total cross section of (6'). 

Thus we conclude that it is reasonable to use the OPE model for calcu- 
lation of the high-energy DCX cross section. 

b) Tests of the assumptions 

Now it is interesting to use the OPE model calculations for checking of 
the simplifying assumptions of our approach to DCX formulated in Ref. [4] and 
discussed above. 

The real and imaginary parts of forward (z = 1) and backward (z = — 1) 
invariant amplitude, A JT + n -(M 2 , z), were calculated using Eqs.(12) - (14). As 



(for t 



max — 
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can be seen from Fig. 11, ReA n + 7T -(M 2 , z = —1) oscillates but it is not small 
compared to ImA 7T + 7r -(M 2 , z = —1) so its contribution should be taken into ac- 
count at energies ~ 1 GeV . This means that an integral over Re^^ (M 2 , z = 
— 1) can not be neglected in this energy range and thus the same is true for 
integral over the large semicircle. The backward parts of the amplitude ap- 
peared to be rather different from the forward ones within various intervals of 
M 2 . The actual ranges of M 2 for a set of initial energies for experimental data 
avaliable on the reaction (6') are given in the Table 2. 

In the Fig. 12 we present the energy dependence of the cross section 
dcrJxjx/dQ, calculated according Eq.(7) at M^ ax \uQ for the following cases: 

(1) ImA'l ™^ is obtained through the Eq.(8) for the forward [z = 1) ampli- 
tudes; 

(2) ImA^^T^ i s replaced by \ m A h £ c ^ ard which is obtained according to Eq.(8) 
for the imaginary part of the backward [z = —1) jttt amplitude; 

(3) Same as in (2), but for ReA b n ™ k 2 ™ ard ; 

(4) both imaginary and real parts of the backward amplitude are taken into 
account. 

The above analysis shows that for 7r + 7r~-scattering amplitudes there are 
substantial differencies between integrals of forward and backward scattering 
amplitudes and they can be neglected only at rather low energies, where s-wave 
production dominates. We think however that OPE model overestimates this 
difference, especially for Tq < 1 GeV, as it predicts too low s-wave production 
at these energies (see Fig. 7). So the OPE-result for DCX can be considered as 
a lower bound for this cross section. 



5 Results 

Let us now calculate the cross section for the reaction 



tT + ib O -> 7T+ + X (16) 

in the kinematical range < AT < 140 MeV of the experiment [6] as a sum 
of the cross sections with 7T° and 2jt in intermediate states (see Eq.(3) ), 



_0 



f doDcx ,do- pcx \ , i dcr 1 pcx \ /-,7\ 

( ^r )l4 ° - ( ^r )l4 ° + { ^r )uo - (17) 



Here (da D ° cx /dQ)uo is the cross section of the reaction (16) within the SSCX 
mechanism (FM effect is taken into account) which is shown in Fig. 13 (see 
dashed curve). (da D n cx /dQ)uo is the contribution of the intermediate tttt state 
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for the given kinematical range of AT. For \mA{ZZ according to Eqs.(3)- 
(5) and normalizing at T Q = 0.6 GeV we obtain the dotted curve in Fig. 13. 
The solid curve corresponds to an account of the full backward scattering 
amplitude. It can be considered as a lower bound for the DCX cross section, 
while the dotted curve can be considered as an upper bound for the cross 
section. Experimental point at 1.1 GeV is close to this bound. This is related 
in our opinion to an underestimate of the s-wave contribution in the OPE 
model at these energies noted above. So at higher Tq we expect future DCX 
measurements to be closer to the solid curve. Note, that both curves at T Q < 
1 GeV, where OPE model strongly deviates from experimental data (see Fig. 6), 
are calculated using not theoretical but experimental values (see the first two 
crosses in Fig. 6). 

The contribution from three-pion intermediate state to Eq.(17) estimated 
using the cross section of 7T~p — > urn appeared to be not more than 10% even 
at the highest energy considered. 

6 Comparison with other approaches 

Let us discuss a relation between inelastic rescatterings considered in this 
paper, and other approaches which have been proposed for the DCX process. 

Relation of the present approach to a model of meson exchanged currents 
(MEC) [18] has been discussed in Ref.[4]. The main difference is that in the 
OPE model we treat the 7r~7r + -scattering amplitude as a function of M 2 and 
integrate over this variable while in the MEC model the 7r~7r + -amplitude is 
approximated as a point- like interaction taken at the threshold (M = 2m n ). 
Due to soft pion theorems this amplitude is small and thus leads to small 
modifications of SSCX predictions [7]. We take into account in the OPE model 
both real and imaginary parts of 7r~7r + amplitude in the regions of M 2 where 
they are not small, so it is not surprising that our results for DCX cross sections 
are substantially higher than SSCX predictions (even for a solid curve in Fig. 
13). 

Another possible mechanism of DCX [19] is a production of an intermedi- 
ate A isobar on one proton of a nucleus, 7T~p — > 7r + A~, with its absorption [20] 
on a second proton, A~p — > nn. A characteristic feature of this mechanism is 
a fast decrease of DCX cross section at To > 1 GeV due to energy dependence 
of the a(7T~p — > 7r + A~) as is shown in Fig. 2 by the curve. 

Our study shows that at least up to energies ~ 4 GeV the region of Mjj ~ 
M^ ax is very essential in integrals over M%, and integrals of both imaginary and 
real parts of amplitudes are important. This means that average distances, d, 
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between participating nucleons are relatively small: d ~ ^m^AT) -1 / 2 . These 
nucleons are closely correlated and can be even considered as a 6-quark system 
(as, for example, in Ref.[21]). 

7 Experimental outlook and conclusions 

We have extended the conventional Glauber mechanism for DCX process 
with an account of inelastic rescatterings. In this paper we investigated the 
problems which arise due to a non-diagonal nature of the process using the 
OPE model. This model gives a reasonable description of the reaction n~p — > 
7r~7r + n at energies above 2 GeV. So we think that our predictions for DCX 
process based on this model, which are substantially higher than SSCX results, 
are reliable in this energy region. DCX is not yet studied experimentally at 
energies above Tq = 1.1 GeV. If strong deviations from our predictions will 
be found it will have important implications: either other mechanisms (e.g. 
those mentioned in the previous Section) are important or the phase-shift 
analysis of ttjt scattering [12], which we have used, is not reliable in the large- 
mass region. An interesting information on these problems can be obtained by 
varying experimental interval of AT. Note that for small AT (and in particular 
in the exclusive limit) an integration region in M 2 is small, and results should 
be closer to SSCX predictions. 

Thus experimental study of the energy dependence of DCX in the region 
of Tq =1 -=- 5 GeV can give an important information on dynamics of this 
process. Cross section in this region according to our estimate is not too small, 
and experiments seems feasible. However, at higher energies the possibility to 
measure the cross section in the region of AT = -r- 140 MeV will be limited 
by the energy resolution of the detector. 

In conclusion, it is important to stress once more that the investigation 
of the pion DCX reactions in the intermediate energy region can provide a 
test on a validity of the Glauber-Gribov approach in the case when effects of 
inelastic rescatterings are large. Another class of processes where such inelas- 
tic contributions are very important is heavy ions interactions at very high 
energies. It is known (see, for example, [22]) that these contributions strongly 
modify a usual Glauber picture and an equilibration in the system can appear 
only as a result of such rescatterings effects. 
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43.6 


39.7 


dvDcx/dQ, fib/sr 
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59.6 ±7.4 


43.3 ±5.5 


26.6 ±8.9 



Table 1. 
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M 2 max {s) 


M 2 max {s)\uQ 


0.61 


0.36 


0.36 


0.88 


0.57 


0.51 


1.02 


0.69 


0.60 


1.20 


0.84 


0.69 


1.45 


1.08 


0.84 


1.76 
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1.02 


2.12 


2.02 


1.26 


2.86 


2.5 


1.68 


3.86 


3.2 


2.2 



Table 2. 
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Figure 1: Diagrams contributing to pion double charge exchange on nuclei: (a) sequential 
single charge exchanges (SSCX), i.e. standard mechanism (elastic rescattering) , (b) quasielastic 
rescatterings, (c) inelastic rescatterings. 
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Figure 2: Experimental total cross section of the reaction n p — > Hn and 7T p — > 7r + A 
compilation [5]). 
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Figure 3: Pion double charge exchange on nucleus (the most general diagram). 
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Figure 6: Total cross section of the reaction — > 7i + 7i~n. Stars and circles are experimental 
[5] and OPE cross sections, crosses and triangles are their corresponding truncated values for 
the AT < 140 MeV; empty star is for limited t (see text). 
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Figure 7: Mass-squared distribution of n + n system in the reaction n p — > n + n n [10 
p = 1.343 GeV/c. The solid curve is our calculation. 




Figure 8: Mass-squared distribution of 7r + 7r~ system in the reaction — > n + n~n [14] at 
p = 1.59 GeV/c, open circles correspond to full range of t, open squares are for limited t (see 
text), the solid and dashed curves are our calculation. 
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Figure 9: Mass distribution of 7v + n system in the reaction tt p — > n + n n [15] at 
p = 2.26 GeV/c. The solid curve is our calculation. 
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Figure 10: Mass distribution of 7r + 7r system in the reaction n p — > 7r + 7r n [16] at 
p = 4 GeV/ c. The solid curve is our calculation. 
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Figure 11: Invariant n + n scattering amplitudes within OPE model. 




Figure 12: (daj? cx / 'dQ)uo calculated using ImA(O) in Eq.(7) (empty triangles), lmA(z = — 1) 
(full triangles), ReA(z = —1) (empty circles), and both ReA(z = —1) and lm(z = —1) (full 
squares) . 
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Figure 13: <Jdcx integrated over AT range from to 140 GeV. 



